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Abstract. - The recently discovered superconductivity in MgB2 has captured world attention 
due to its simple crystal structure and relatively high superconducting transition temperature 
Tc = 39K. It appears to be generally accepted that it is phonon-mediated s-wave BCS like 
superconductivity. Surprisingly, the strongly temperature dependent anisotropy of the upper 
critical field, observed experimentally in magnesium diboride single crystals, is still lacking a 
consistent theoretical explanation. We propose a simple single-gap anisotropic s-wave order 
parameter in order to compare its implications with the predictions of a multi-gap isotropic 
s-wave model. The quasiparticle density of states, thermodynamic properties, NMR spin- 
lattice relaxation rate, optical conductivity, and Hc2 anisotropy have been analyzed within this 
anisotropic s-wave model. We show that the present model can capture many aspects of the 
unusual superconducting properties of the MgB2 compound, though more experimental data 
appear to be necessary from single crystal MgB2. 



Introduction. - The moderately high temperature superconductivity with Tc ~ AOK, 
discovered about one year ago in MgB2 ||^, || has stimulated an intense research activity all 
over the world. Superconductivity in this binary compound appears to be due to an electron- 
phonon interaction and compatible with BCS s-wave superconductivity. Thermodynamics 
and contact tunnehng data as well as some theoretical studies indicate that superconductivity 
in MgB2 is one of the rare examples of two-band superconductivity with two energy gaps, 
attached to different sheets of the Fermi surface [|-10|. On the other hand, the two gap 



model appears not to be able to cope with a strongly anisotropic upper critical field in c- 



axis oriented MgB2 films and more recently in single crystals of MgB2 11-16|. Indeed, the 
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Fig. 1 - Anisotropic s-wave order parameter Eq. (|l|) for a — 10 
Fig. 2 - Density of states for a— 5, 10, and 20. 

anisotropic s-wave model can describe in principle an anisotropy of the upper critical field and 
an angular dependence of i?c2 including deviations from the Ginzburg-Landau prediction with 
anisotropic mass term Indeed, some of the STM data appear to be more consistent 

with an anisotropic s-wave model | [T9| ]. However, more telling is the temperature dependence 
of the ratio of i?c2(^) ( ^ = T/Tc) and i?c2(0 where the superfix a and c denotes a field 
oriented parallel to the a-axis or the c-axis, respectively. The ratio ^{t){= "2/^02) increases 
as temperature decreases as has been observed in Refs. [p"2|-p^. This clearly indicates that 
A(k) has to have oblate rather than the prolate form suggested earlier |jl^,|l8|. Otherwise 
7(i) would decrease with decreasing temperature. 

The object of this paper is to propose an oblate A(k) gap anisotropy and to see whether 
this order parameter can describe the two gap features observed experimentally. Indeed 
choosing one adjustable parameter, which determines the ratio of Amin/Amax, we can describe 
reasonably well the density of states measured by STM and the specific heat data of Wang 
et al [|| . Also some peculiarities of the NMR relaxation rate and optical conductivity can be 
described. Then we shall go on to study the upper critical field for H || c and H || a. 

As we shall see later, we can describe i?c2(0 reasonably well by choosing our parameter 
a > 5 (corresponding to a ratio of A,„in/Amax < 0.4, see below). On the other hand, the strong 
increase of -^^2(0 appears to be somewhat difficult to fit consistently with other experimental 
data within our model. Of course in the real single crystals the electronic mean free path is 
rather short / ~ 60 - 70A |l|-|l|,|2^,|2|. Therefore in a more realistic analysis it is necessary 
to consider impurity scattering as well. 

Thermodynamics. - We consider an anisotropic BCS model for superconductivity in 
MgB2 with an order parameter given by 



A(k) = A 



1 



(1) 



where the parameter a determines the anisotropy, z = cos{9) and 9 is the polar angle with 
respect to the c-axis. In Fig. Hthe anisotropic s-wave order parameter is plotted in momentum 
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Table I - Gap ratios and specific heat jumps for a— 5, 10, 20, and 40. 



a 






AC 


5 


2.29 


0.93 


1.11 


10 


2.46 


0.74 


0.93 


20 


2.61 


0.57 


0.74 


40 


2.78 


0.43 


0.57 



space for a ~ 10 displaying a pancake-like shape. Its maximum value A,nax = A lies in the 
a6-plane, while in c-axis direction it assumes its minimum value of Amin = A/vT+a. In the 
following we will explore the consequences of anisotropic s-wave superconductivity in MgB2 
within the framework of the weak-coupling BCS theory. 

We have solved numerically the weak-coupling gap equation for the state Eq. ffl) 



/>OC ( 

Jo '^'W'^'^y 




tanh ^ ^ ^ - -tanh ^ = (2) 



Here, f{z) = l/\/l + az^ and (• • •) = dz - ■ ■ denotes an angular average over the variable 
z. In Table | we tabulated gap ratios found from the solution of Eq. (g) for a= 5, 10, 20, 
and 40. With increasing value of a an increasing gap ratio of the maximum gap Amax/^c is 
found, while the minimum gap ratio appears to decrease. 

In Fig. H we show the corresponding density of states for these values of a calculated from 

The data in p9| resembles our density of states if we choose a ~ 3. Further the peak in 



the density of states in |19| gives a maximum gap ~ 8 meV consistent with our model. Also 
the authors [|l^ deduced a momentum dependence of A (k), which is topologically the same as 
the earlier model [^,|l^. On the other hand, the data in |2^ exhibits a two gap like structure 
very different from Fig. ^j. Further, the maximum gap 7 meV is somewhat smaller then the 
one expected from the present model. However, the double gap like feature is inherent to the 
anisotropic gap, since the two extrema associated with the single gap are always visible jr7| , p^ . 

In Fig. ^ the specific heat is shown for values of a= 10, 20, and 40 along with the data 
of Wang et al. [^|^. The specific heat jumps found from our model are also tabulated in 
Table ||. They decrease with increasing a. The closest agreement with the experimental data 
is obtained for a — 20. However, fitting of the structures seen in the experimental data would 
require some fine tuning of the angular dependence in Eq. (^. 

In Fig. ^ we show the normalized NMR spin-lattice relaxation rate obtained from our 
state. With increasing anisotropy a the Hebel-Slichter peak is reduced. Roughly, the size of 
the peak found here is consistent with recent NMR data on polycrystalline MgB2 by Kotegawa 
et al. (2|. 

It has recently been noted that the gap seen in optical conductivity appears to be very 
small I2J, We calculated the anisotropic analog of the Mattis-Bardeen conductivity 
averaged over all directions, which should correspond to the response of a polycrystalline thin 
film sample as has been used in Refs. [^^. For our anisotropic state the conductivity a{uj) 



is obtained using the following two functions: 
ii[x) = I az I ay- 



'0 
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hix, Zq) 



dz / dy—= 

zo Jf{z) v(y 



P{z)-y{x~y) 



P{z))[{y-xY-P{z)] 



: tanh 



2T J 



Using these the optical conductivity reads 



4A2 



for W < 2Amin 

for 2A,nin < LO < 2A 

max 

for UJ > 2Amax 



(4) 



The result for a = 20 and different reduced temperatures is shown in Fig. ||. Indeed it appears 
that the optical conductivity is mainly dominated by the minimum gap within our model, 
because the absorption threshold starts at 2Amin- Thus, an anisotropic gap scenario can 
also account for the small gap seen in the optical conductivity. Most useful would be optical 
conductivity studies on single crystals which could resolve the c-axis and a6-plane response in 
order to see more clearly, what the anisotropy looks like. 

Upper critical field. - Recent experiments on the upper critical field in MgB2 single crys- 
tals have reported a strong temperature dependent anisotropy of Hc2 and an unusual upward 
curvature of the critical field parallel to the ab plane. An anisotropy of the superconducting 
pairing can in principle give such an effect. We investigate the temperature dependence of 
both Hc2, parallel and perpendicular to the crystal c-axis of MgB2. 

The general equation for the upper critical field is derived from the gap equation. For 
unconventional superconductors it can be treated variationally, some technical details can be 
found elsewhere ||27| . 
(a) H II c ^ 

The equation for the upper critical field parallel to the c-axis is given by 



\nt = 



du 
sinh I 



l-{f)-'{eM-pu'{l-z'))f)\, 



(5) 
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Fig. 5 - Mattis-Bardeen conductivity as a function of normalized frequency a;/ A for a = 20 and 
reduced temperatures t — T/Tc of 0.9, 0.8, 0.6, 0.4, and 0.2 (from top to bottom), respectively. 

Fig. 6 - Normalized upper critical field hc2 as a function of reduced temperature t for a = 10 and 
o = 100. The upper two curves are for field direction within the a6-plane, while the lower two curves 
are for field in c-axis direction. For comparison also the isotropic result (a — 0) is shown. The data 
points are taken from Ref. ||l^, where we have normalized them to the slope at Tc, taking Tc = 37. 3K, 
'^"dx'^' \Ta - -0.435T/K, and ^^f^k, =i -0.124T/K. 



where p = v'^eHc2{t) /2{2ttT)'^ , and Va is the Fermi velocity within the a5-plane. 

In Fig. 1^ the numerical solution of this equation, normalized by its derivative at T^, 
/ic2(T) = H^2{T)/{-T^dHc2iT)/dT)\T,, is plotted along with the /ie2(r)-curves for the ab- 
direction, the corresponding curve for the isotropic case, and normalized experimental data by 
Angst et al. |l3j. In the limit T ^ for a = 10 we have /ic2(0) — 0.63349 which is somewhat 
smaller than the corresponding value in the isotropic s-wave superconductor /ic2(0) — 0.72726. 
In the vicinity of Tc the upper critical field exhibits naturally a rather linear temperature de- 
pendence within the weak-coupling BCS theory. This behavior is observed for both anisotropic 
s-wave and conventional superconductors. Note, that our results for the upper critical field 
parallel to the c-axis for a > 5 are in good correspondence with the experimental data jist . 
The best fit to these data is obtained for a w 10. However, we cannot reliably extract a value 
of a from these data, because h^2i^) for a — + oo is saturating at a value of 0.59054 and the 
data are already close to this limit. Thus, also even higher values of a might still be consistent 
with the c-axis data. 

(b) H II a 

The derivation of the equation for the temperature dependent upper critical field in the 
plane is more involved. The problem is that a mixing of higher Landau levels takes place p^ . 
We choose a variational wavefunction suggested in Ref. |^, Eq. (36), which corresponds to 
a distorted Abrikosov ground state. Due to the anisotropy of our state such a distortion is 
expected for the field direction perpendicular to the c-axis. The distortion can be varied using 
a parameter a, which has to be determined by a variational principle. Using this method 
we arrive at the following equation for the temperature dependence of the upper critical field 
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Pa{t) = VaVceHc2it)/2{27rTf 

- ^ i" ssl^f - <^'>" (»■> (i" - ™^ * + "'=0 } }• 

here (•■■) = Jq dz J^^ 5^ ' ' ' ^.nd is the Fermi velocity in c-axis direction. The parameter 
a has to be determined from a minimization of the right-hand side of Eq. (H), taking its 
derivative to be zero. 

We have compared our results with a second, independent calculation based on a Landau 
level expension up to the fourth Landau level following the method of Luk'yanchuk and 
Mineev . For the parameter range a < 20 both methods were in close agreement with each 
other. However, for a ^ 20 the Landau level expansion started to give significantly smaller 
values for i?°2i indicating that inclusion of even higher Landau levels becomes necessary, and 
showing that the distorted Abrikosov wavefunction is a better variational solution in this case. 

The numerical solution of Eq. (|^) obtained from a numerical optimization of a is shown in 
Fig. H In the zero-temperature limit we obtain = 0.87309 for a = 10 and = 1-18025 
for a — 100. The curve for a = 100 possesses a slight upward curvature in agreement with 
the experimental data. From our comparison with the experimental data [ p^ we extracted 
the values -dH^2(T)/dT ~ 0.435 T/K and ^dH^2{T)/dT ~ 0.124 T/K. Choosing a = 10, 
from these we deduce average Fermi velocities Va — 2.7 • 10*" cm/sec, Vc — 1.3 • 10'' cm/sec and 
the ratio Vc/va — 0.48. The last ratio is very consistent with what was obtained earlier |]l^ 
using the data in |0 . Also the absolute values of these Fermi velocities Va and Vc are roughly 
consistent with bandstructure calculation results provided to us by A. Yaresko. 

This rather large value of a = 100 appears to give a reasonable fit to the a6-plane data 
of Hc2- But unfortunately, this corresponds to a ratio of Amin/Amax — 0.1, which appears 
to be inconsistent with the values we discussed above. Here we can just speculate about a 
possible resolution of this problem: in these single crystals the mean free path is relatively 
short, being about the order of the coherence length j2l],^, which puts these systems in an 
intermediate region between the clean and the dirty limit. Therefore a study of the influence 
of impurity scattering would be necessary. It is possible that the upper critical field anisotropy 
is much less sensitive to impurity scattering than the energy gap and the density of states, for 
example. Such an effect could possibly resolve this problem and is left for a future work. 

Conclusions. - We have described a model of anisotropic s-wave superconductivity and 
compared our theoretical predictions with the experimental data of MgB2. We show that our 
model can capture many aspects of the two gap model and is consistent with the experimental 
data of single crystals MgB2. Also we stress that so far our model is the only one, that can 
describe the strongly temperature dependent anisotropy in the upper critical field. Unfortu- 
nately, strongly different values of our anisotropy parameter a had to be used for different 
experimental properties. On the other hand, we have ignored the effect of impurity scattering 
in the present analysis. In particular the mean free path / ~ 60 — 70A reported for single 
crystals of MgB2 implies that we are in the intermediate regime (neither in the clean limit 
nor in the dirty limit). Therefore in a more realistic analysis, it will be very important to 
incorporate the effect of impurity scattering. Clearly, we also need more precise measurements 
on single crystal MgB2. 

* * * 
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